Hydrogen-enhanced local plasticity in aluminum: an ab initio study 
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Dislocation core properties of Al with and without H impurities are studied using the Peierls- 
Nabarro model with parameters determined by ab initio calculations. We find that H not only 
facilitates dislocation emission from the crack tip but also enhances dislocation mobility dramat- 
ically, leading to macroscopically softening and thinning of the material ahead of the crack tip. 
We observe strong binding between H and dislocation cores, with the binding energy depending 
on dislocation character. This dependence can directly affect the mechanical properties of Al by 
inhibiting dislocation cross-slip and developing slip planarity. 
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The mechanical properties of solids are not only a func- 
tion of their intrinsic atomic structure, but also of the 
environment in which they exist. Small amounts of im- 
purity atoms in the environment can drastically change 
the response of a solid to external loading. Hydrogen em- 
brittlement is one of the most important and well-studied 
processes of environmental degradation of materials. Af- 
ter over a century of studies, the definitive mechanism 
of H embrittlement still remains unknown although sig- 
nificant progress has been made toward developing a 
detailed understanding of the problem. Three general 
mechanisms for H embrittement have been proposed: (1) 
Formation of a hydride phase; (2) H-enhanced local plas- 
ticity; and (3) H-induced grain boundary weakening |Q. 
The underlying atomic processes and the relative impor- 
tance of the three mechanisms, however, remain uncer- 
tain and controversial. In this paper, we study a key as- 
pect of H embrittlement, referred to as H-enhanced local 
plasticity (HELP) ||, which has attracted much atten- 
tion recently. There seems to be overwhelming evidence 
that plasticity is a fundamental contributor to H embrit- 
tlement of virtually all susceptible metals and alloys, in- 
cluding Al, the model material we consider in this paper 
KLpH)- The most dramatic experimental observation is 
that H increases dislocation mobility significantly under 
constant stress. This H-enhanced mobility is observed 
for screw, edge and mixed dislocations as well as for iso- 
lated dislocations and dislocation tangles pM]. Other 
important experimental results include the observations 
of slip planarity and strong binding between H and 
dislocation cores 

In contrast to the vast body of experimental evidence, 
theoretical studies of the HELP mechanism are scarce. 
To our knowledge, the only theoretical investigation of 
the HELP mechanism is based on finite-element modeling 



of elastic interactions among dislocations in the presence 
of H H . Although this approach was able to explain the 
H-enhanced dislocation mobility through the H-induced 
"elastic shield" effect, it completely ignores the effects of 
H on dislocation core structure and contains no informa- 
tion for the relevant atomic processes, which are critical 
for a complete understanding of the phenomenon. The 
purpose of our work is to provide a comprehensive un- 
derstanding of the aforementioned experimental results 
based on an ab initio theoretical framework, which hope- 
fully can set the stage for an atomistic theory of H em- 
brittlement. 

In this paper, we have used the recently developed 
Semidiscrete Variational Peierls-Nabarro (SVPN) model 
|^| |To| with ab initio determined 7-surfaces jll| and elas- 
tic constants. This approach has been shown to predict 
reliable dislocation core properties by comparing its pre- 
dictions to direct atomistic simulations based on the same 
force law as that used for the 7-surface calculations |||| . 
The uniqueness of the approach when combined with ab 
initio calculations for the energetics is that it produces 
essentially an atomistic simulation for dislocation core 
properties without suffering from the uncertainties asso- 
ciated with empirical potentials. For example, atomistic 
simulations based on EAM potentials predict that dis- 
locations in Al will dissociate into partials [fllLOJ, while 
experimentally no such dissociation is observed p4[ , a 
discrepancy which can be traced to the fact the EAM 
potentials underestimate the intrinsic stacking fault en- 
ergy [§• 

In the SVPN approach, the equilibrium structure of 
a dislocation is obtained by minimizing the dislocation 
energy functional 

Udisl = Udastic + U m i s fit + U st ress + Kb hlL, (1) 
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with respect to the dislocation displacement density 
|p|-|To|. The 7-surface determined from ab initio calcula- 
tions enters the equation through the U m i S fit term. We 
identify the dislocation configuration-dependent part of 
the elastic energy and the misfit energy as the core en- 
ergy, U core = V r e i as tic + U mis f it , while the configuration- 
independent part of the elastic energy Kb 2 h\L is excluded 
because it is irrelevant in the variational procedure and it 
has no contribution to the dislocation core structure. The 
response of a dislocation to an applied stress is achieved 
by the minimization of the energy functional at the given 
value of the applied stress. An instability is reached when 
an optimal solution for the density distribution no longer 
exists, which is manifested numerically by the failure of 
the minimization procedure to convergence. The Peierls 
stress is then defined as the critical value of the applied 
stress giving rise to this instability. 

In order to study the interaction between H and dis- 
locations, we have carried out ab initio calculations of 
the 7-surface with H placed at interstitial sites of the Al 
lattice. Specifically, we have used a large supercell con- 
taining six Al layers in the [111] direction with four atoms 
per layer to simulate a low concentration of H impurity 
(~ 4 at.%). Two sets of calculations were performed with 
H placed at the octahedral and the tetrahedral sites in 
Al. The ab initio calculations are based on the pseu- 
dopotential |l5| plane-wave method with the generalized 
gradient approximation to the exchange-correlation func- 
tional |l(|. A kinetic energy cutoff of 12 Ry for the plane- 
wave basis is used and a fc-point grid consisting of (8,8,4) 
divisions along the reciprocal lattice vectors is sampled 
for the Brillouin zone integration. In our calculations, 
the H atom is allowed to relax freely according to the 
Hellmann-Feynman force acting on it while Al atoms are 
restrained to move only along the [111] direction. Vol- 
ume relaxation is performed for each sliding distance to 
minimize the tensile stress on the supercell. Table I lists 
the various stacking fault energies obtained for the pure 
Al and the Al+H systems. Although the octahedral site 
is more stable for H without sliding, the energy for the 
tetrahedral site can be lower for certain sliding vectors. 
Therefore the energies in Table I and Fig. 1 are given 
with respect to the lowest values of the two sites for each 
sliding. From Table I we find that H decreases the energy 
of the various stacking faults by up to 50%, which could 
change dislocation properties significantly. For example, 
one can predict that dislocations would be emitted more 
easily from a crack tip when H is present based on the 
lowering of the unstable stacking fault energy along the 
[121] direction (l7). In fact, this prediction is in accor- 
dance with atomistic simulations for Ni, where disloca- 
tions are found to be emitted more rapidly from a crack 
tip in the presence of H |l^] . 

The 7-surface for both the pure Al and the Al+H sys- 
tems is shown in Fig. 1(a) and (b) respectively. The 
energy along various high symmetry directions is deter- 



mined from ab initio calculations while the rest of the 7- 
surface is obtained by fitting a symmetrized polynomial 
basis, which guarantees the correct rotational symmetry 
of the 7-surface on the (111) plane. Comparing the two 
7-surfaces, we find an overall reduction in energy in the 
presence of H. In order to understand how H affects the 
elastic properties of Al, we also calculate the various elas- 
tic constants for both the pure Al and the Al+H systems 
using ab initio approaches. We used a smaller unit cell 
with only four Al atoms and one H impurity to facilitate 
these calculations, and find that, even at such a high con- 
centration of H, the prelogarithmic elastic energy factors 
entering the model J9| change only by a few percent at 
most |18|| . 

Having determined all the necessary parameters that 
enter the model, we can start to investigate the core prop- 
erties for various dislocations in Al. We have studied 
screw (0°), 30°, 60° and edge (90°) dislocations. The 
Peierls stress for these dislocations is shown in Table II. 
We find that the Peierls stress, which is the minimum 
stress to move a stationary dislocation, is reduced by 
more than an order of magnitude in the presence of H. 
Moreover this H-enhanced mobility is observed for screw, 
edge and mixed dislocations, which is compatible with 
the experimental results. Not only is our calculation the 
first successful attempt to explain HELP associated with 
individual dislocation behavior, it also invalidates the 
perception that the only plausible explanation of HELP 
must be based on elastic interactions among dislocations 
and that dislocation-lattice interaction is not important 
[jij. On the contrary, we find the dislocation-lattice in- 
teraction to be very important and that the dislocation 
core structure is responsible for the observed H-enhanced 
dislocation mobility. Our result, however, does not ex- 
clude the possibility that dislocation-dislocation interac- 
tions do play a role in HELP. 

Other important experimental results that any the- 
ory of H cmbrittlcmcnt must explain are the observed 
H trapping in dislocation cores and the H-induced slip 
planarity. To this end, we have calculated the energetics 
of dislocations for the pure Al and the Al+H systems. 
The results are shown in Table II. We remind the reader 
that E core listed here is only the displacement density 
dependent part of the total energy, excluding the large 
positive elastic term Kb 2 \nL which renders the total en- 
ergy positive j^] . The binding energy of H is defined as 
the difference between dislocation core energies with and 
without the presence of H. From these results we find that 
there is strong binding between H and dislocation cores 
and that H is attracted (trapped) to dislocation cores to 
lower the core energies. More importantly, we find that 
the binding energy is a function of dislocation character, 
with edge dislocations having the greatest and screw dis- 
locations having the lowest binding energy, respectively. 
For a mixed dislocation, we find that the binding energy 
increases with the amount of edge component. These re- 
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suits have significant consequences for the experimentally 
observed inhibition of cross-slip and slip planarity in the 
presence of H. Experiments show that in the presence of 
H slip is confined to the primary slip plane and disloca- 
tion cross-slip is restricted, which leads to the observed 
slip planarity. On the other hand, removal of H from the 
sample allows a dislocation segment on the primary plane 
to reorient towards the screw orientation and to complete 
the cross-slip precess j2(J. It was believed from the ex- 
perimental evidence that any process, such as cross-slip, 
that increases the screw component length at the expense 
of the edge component will be made more difficult by the 
presence of H |pof . Our results for the binding energy 
of H in dislocation cores are the key for understanding 
these experimental observations: Since the edge disloca- 
tion has twice the binding energy as the screw, it will 
cost much more energy for an edge dislocation in the 
presence of H to transform to a screw dislocation in or- 
der to cross-slip. In the same vein, it costs more energy 
for a mixed dislocation to transfer its edge component 
to a screw dislocation for cross-slip. In other words, H 
can stabilize the edge component of a mixed dislocation 
and inhibit the dislocation cross-slip process, leading to 
the observed slip planarity. It is also interesting to point 
out that our results for the binding energy of H are in 
qualitative agreement with those from atomistic simula- 
tions of H in Ni, where the binding energies of H to edge 
dislocations and screw dislocations are found to be 0.12 
eV and 0.09 eV per H atom, respectively |2l| . 

To further understand these results, we calculate the 
displacement density for the various dislocations which is 
presented in Fig. 2. We find that although the intrinsic 
stacking fault energy is reduced by almost 40% in the 
presence of H, dislocations are not dissociated into par- 
tials [^U . This result is important because it discredits 
the argument that the observed restriction of cross-slip 
may result from the increased separation of partials due 
to the reduction of the stacking fault energy (23). In 
fact, our result is consistent with the experimental find- 
ings that no widely spaced dislocation partials were ever 
observed in TEM studies of Al containing large concen- 
trations of H pOR . Therefore, we believe that the ob- 
served slip planarity and restriction of cross-slip is due 
to the strong binding between H and the edge compo- 
nent of dislocations, not only for Al, but possibly also for 
other metals ||||§(J. 

Based on the experimental observations and the re- 
sults we obtained in this study, we identify the following 
processes which we believe could lead to the H-induced 
fracture: (a) H is attracted to crack tips: Application 
of external stress produces local concentration of tensile 
stress in the vicinity of cracks, which attracts H since H 
prefers to stay in slightly enlarged interstitial sites, (b) 
Thinning and softening: These segregated H impurities 
at the crack tip can facilitate dislocation generation and 
enhance dislocation mobility, which will lead to extensive 



plastic deformation in front of the crack, causing thinning 
and softening of the material ahead the crack, (c) Lower- 
ing fracture stress: The thinning and softening processes, 
along with the localization of slip due to the inhibition 
of cross-slip, allow the crack to propagate at lower stress 
levels, prior to general yielding away from the crack tip. 
Although our calculations concern Al, we believe that 
the results are also applicable to other metals. Further- 
more we should point out that the above processes are 
not likely to be the only mechanism operating in H em- 
brittlement j^] , but are the most relevant processes based 
on HELP theory. 

In summary, we have studied the interaction of inter- 
stitial H with dislocations in Al using the SVPN model 
with ab initio determined parameters. We find that H 
can lower the Peierls stress of dislocations by more than 
an order of magnitude, strongly supporting the HELP 
model. We also find there is a strong binding of H to 
dislocation cores and that the binding energy is a func- 
tion of dislocation character. Our findings explain the 
experimentally observed restriction of the cross-slip pro- 
cess and slip planarity in the presence of H. The dislo- 
cations do not dissociate into partials even though the 
intrinsic stacking fault energy is reduced by 40% in the 
presence of H. Finally, based on the results of our calcu- 
lations, we identify the processes which could lead to H 
cmbrittlement. 
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TABLE II. Peierls stress (a p , MPa), core energies (E core , 
eV/A) for the four dislocations in the pure Al and the Al+H 
systems and binding energy (Et, eV/atom) for the four dis- 
locations. 







screw 


30° 


60° 


edge 




Al 


254 


51 


97 


3 




Al+H 


1.7 


1.2 


1.4 


0.3 


E 


Al 


-0.08 


-0.11 


-0.17 


-0.20 




Al+H 


-0.14 


-0.18 


-0.27 


-0.32 


E b 




0.06 


0.07 


0.10 


0.12 



FIG. 1. The 7-surfaces (J/m 2 ) for displacements along a 
(111) plane for (a) pure Al and (b) Al+H systems. The cor- 
ners of the plane and its center correspond to identical equi- 
librium configurations, i.e., the ideal lattice. The two energy 
surfaces are displayed in exactly the same perspective and on 
the same energy scale to facilitate comparison of important 
features. 



FIG. 2. Dislocation displacement density for four disloca- 
tions (clockwise): screw (0°), 30°, 60° and edge (90°) for the 
pure Al (solid lines) and the Al+H (dashed lines) systems. 



TABLE I. Fault vectors and energies (J/m 2 ) for four im- 
portant stacking faults of the pure Al and the Al+H systems. 





Vector 


Al 


Al+H 


Intrinsic stacking 


1/6[121] 


0.164 


0.103 


Unstable stacking 


1/10[121] 


0.224 


0.126 


Unstable stacking 


1/4[101] 


0.250 


0.132 


Run-on stacking 


1/3[121] 


0.400 


0.310 
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